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(57) Abstract 

A device (10) that detects and reports the presence of an 
acute cardiac ischemic condition in a patient includes sensing 
electrodes (12, 14, 16, 18, 20, 22, 24, 26, 28, and 30) placed 
on a partient for acquiring multiple (e.g., twelve) lead ECG data 
from the patient. TTie device evaluates the acquired ECG data 
by analyzing global features derived from the ECG data. Global 
features include projection coefficients calculated from projecting 
a concatenated vector of representative heartbeat data onto sets of 
basis vectors that define an acute cardiac ischemic ECG subspace 
and a non-ischemic ECG subspace. A classifier evaluates the global 
features to determine whether an acute cardiac ischemic condition 
is detected. In a further embodiment of the invention, one or 
more classifiers evaluate local features, such as local morphological 
features and patient clinical information, in addition to the global 
features to determine whether an acute cardial ischemic condition 
is detected. The device (10) may be configured with an adjustable 
sensitivity/specificity operating point. The result of the evaluation 
is reported to the user of the device. 
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METHOD AND APPARATUS FOR DETECTING A CONDITION 
ASSOCIATED WITH ACUTE CARDIAC ISCHEMIA 

Related Applications 
This application claims the benefit of U.S. Provisional Application Serial 
5 No. 60/083,722 filed April 30, 1998, and U.S. Provisional Application Serial 
No. 60/100,391 filed September 15, 1998. 

Field of the Invention 
The present invention relates generally to. the analysis of cardiac electrical 
activity, and more specifically, to the evaluation of electrocardiogram data to detect 
10 and report cardiac abnormalities. 

Background of the Invention 
A variety of physiological processes are electrically mediated and produce 
asso ciated electrical signals. For example, the sinoatrial node in a human heart 
generates an electrical pulse that triggers the remainder of a heartbeat in a normally 
15 functioning heart. This pulse propagates through the heart's normal conduction 
pathways producing electrical signals that can be observed on the surface of a 
patient's body. Monitoring and analysis of such electrical signals have proved 
beneficial in evaluating the function of a patient's heart, including the detection of 
conditions associated with acute cardiac ischemia. 
20 Monitoring of a patient's cardiac electrical activity is conventionally performed 

using a 12-lead electrocardiogram (ECG) system that includes a monitor and ten 
electrodes attached to a patient. A conventional 12-lead ECG system monitors the 
voltages sensed by the ten electrodes and generates twelve combinations of these 
voltages to produce the "leads" required by the 12-lead ECG system. Of the ten 
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electrodes in a 12-lead ECG system, four are "limb" electrodes typically placed on or 
near each of a patient's four limbs, and six are "precordial" electrodes positioned on 
the patient's chest over the heart. As an electrical impulse propagates through the 
heart, the monitor repetitively measures the voltages sensed by the electrodes. 
Although the electrodes collectively monitor the same heartbeats, the electrodes sense 
different voltages due to their placement with respect to the patient's heart. A time 
sequence of monitored voltages is used to produce ECG lead data. An ECG monitor 
typically plots this data to -provide graphical waveforms representing the heart's 
electrical activity for each lead being monitored. 

An example of an ECG waveform is shown in FIGURE 1. For purposes of 
analysis, an ECG waveform produced over a time interval corresponding to one 
cardiac cycle or "heartbeat" is divided into a number of waves. The portion of a 
waveform representing depolarization of the atrial muscle fibers is referred to as the 
"P" wave. Depolarization of the ventricular muscle fibers is collectively represented 
15 by the "Q," "R," and »S" waves of the waveform. The portion of the waveform 
representing repolarization of the ventricular muscle fibers is known as the "T" wave. 
Between heartbeats, an ECG waveform returns to an "isopotential" line. FIGURE 1 
also illustrates selected fiducial points labeled "q", "j", »tl", and "t2". Fiducial points 
define the boundaries of selected features and are used in measuring characteristics of 
20 an ECG waveform, such as the start and end of a heartbeat and the elevation of the 
ST portion of a heartbeat. The "q" point shown in FIGURE 1 represents the start of 
the Q wave, the "j" point represents the end of the QRS complex, the "tl" point 
represents the start of the Twave, and the "t2" point represents the end of the 
T wave. 

25 As noted, an analysis of a patient's ECG may assist in detecting acute cardiac 

ischemia in the patient. As a matter of background, acute cardiac ischemia is a 
condition that arises from chronic or sudden onset of deprivation of blood, and hence 
oxygen, to muscles of the heart. If an ischemic condition is severe or prolonged, it 
can result in irreversible death or damage to myocardial cells (i.e., an infarction). ' A 
30 chronic cardiac ischemic condition, angina, is typically caused by narrowing of the 
coronary arteries due to spasms of the wall muscles or partial blockage by plaques. A 
sudden cardiac ischemic condition may be caused by a clot blocking the passage of 
blood in the coronary arteries. Symptoms of a cardiac ischemic event may include 
chest pain and pain radiating through the extremities, but not all such events present 
these symptoms. Current medical intervention for severe acute ischemic events 



35 



WO 99/55226 PCT/US99/0931 1 

-3- 



includes the administration of a class of drugs called thrombolytics that dissolve clots 
in the occluded coronary artery, and emergent PTCA, a medical procedure that opens 
the artery by inflating a balloon inside the clot to make a passage for circulation. 

The amount of damage done to a heart by an ischemic event depends, in part, 

5 on the amount of time that lapses before treatment is provided. Therefore, ECG data 
should be evaluated as early as possible so that functional changes associated with 
cardiac ischemia can be detected and reported as early as possible. With early 
detection of acute cardiac ischemia, appropriate treatment can take place as early as 
possible and, thereby, maximize the preservation of myocardium. The American 

10 Heart Association recommends that a patient with suspected acute cardiac ischemia 
be evaluated by a physician using a 12-lead ECG within ten minutes of arrival at a 
hospital emergency department. Unfortunately, outside of a hospital, highly trained 
medical personnel are not always available to meet a patient's immediate needs. Quite 
often, a first responding caregiver to a patient in the field may not be competent in 

15 evaluating ECG waveforms to detect acute cardiac ischemic events. A need, 
therefore, exists for a device that not only obtains ECG data, but also quickly 
evaluates and automatically produces a preliminary diagnosis that an acute cardiac 
ischemic condition has been detected. 

Traditionally, acute cardiac ischemic conditions are detected by a physician's 

20 visual evaluation of 12-lead ECG waveforms. A physician typically selects one or 
more leads in the 12-lead ECG and makes an initial assessment by comparing selected 
features of the patient's ECG waveforms to equivalent features of other persons' ECG 
waveforms that are representative of various abnormal conditions. A physician may 
also look at the patient's ECG waveforms over time and evaluate any changes in 

25 waveform shape. A number of waveform features have been identified as usefiil in 
diagnosing acute cardiac ischemic conditions. Customarily, a physician observes the 
extent to which the ST portion of a waveform exceeds the isopotential line (i.e., the 
ST elevation) and uses this information to determine if an acute ischemic event has 
occurred. Nevertheless, because of the subtleties involved in evaluating ECG 

30 waveforms, even highly trained individuals often fail to correctly diagnose an acute 
cardiac ischemic event when evaluating an ECG using traditional features alone. 
More subtle, globally distributed ECG features remain undetected. A need, therefore, 
exists for more accurate ways of detecting and reporting acute cardiac ischemic 
events. 
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In recent years, there have been efforts to develop enhanced ECG waveform 
interpretation based on computer analysis. Conventional computer processes used for 
ECG waveform analysis are based on heuristics derived from the experience of expert 
physicians. Such processes implement rules that attempt to simulate an expert 
physician's reasoning but perform no better than the expert. In practice, many such 
processes perform more poorly than human expert evaluation. 

Furthermore, when a conventional heuristic process is used, it is difficult to 
choose an optimal operating point for the device in terms of sensitivity (i.e , detects 
true positives) and specificity (i.e., avoiding false positives). A device tuned to be 
more sensitive is typically less specific, while a device tuned to be more specific is 
typically less sensitive. A typical sensitivity/specificity tradeoff is illustrated by a 
receiver operating characteristics (ROC) curve, an example of which is shown in 
FIGURE 15. Using a conventional heuristic process, it is difficult to make the 
sensitivity/specificity tradeoff. explicit; thus, the selection of an operating point on the 
ROC curve is often made in a suboptimal, ad hoc manner. As such, there is a need for 
an apparatus and method that can provide better performance in terms of sensitivity 
and specificity and further provide for selecting a sensitivity/specificity tradeoff in a 
more systematic way. 

The present invention provides a method and apparatus having such features, 
20 as well as addressing other shortcomings in the prior art. 

Summary of the Inventinn 
In accordance with the present invention, globally-distributed features in a 
patent's ECG are extracted and analyzed to detect an acute cardiac ischemic 
condition. First, representative heartbeat data for a patient is calculated from the 
patient's raw ECG data for each of the leads obtained. The representative heartbeat 
data from each lead is concatenated to form a vector of representative heartbeat data 
for the set of measured leads. The concatenated vector of representative heartbeat 
data is then mathematically projected onto predetermined basis vectors that define 
signal subspaces of ECGs exhibiting acute cardiac ischemic and non-ischemic 
30 conditions. The resulting projection coefficients are "global features" that are 
evaluated by a classifier to determine whether the global features are indicative of an 
acute cardiac ischemic condition. Analysis of global features according to the present 
invention has been found to improve sensitivity while maintaining high specificity as 
compared to an analysis performed by conventional ECG analysis. 
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In accordance with further aspects of the invention, local features may be 
analyzed in addition to global features to further enhance diagnostic accuracy. Local 
features include measures of morphological features defined for a local time interval 
less than the duration of a representative heartbeat. An example of a local 

5 morphological feature is a measure of ST elevation derived from individual lead data. 
Clinical information, such as a patient's age and sex, may also be included in the 
analysis as local features. The local features and global features are either separately 
evaluated with the results combined, or evaluated in combination to produce a single 
result, yielding a final decision of whether an acute cardiac ischemic condition is 

1 0 present in the patient. 

In accordance with other aspects of this invention, a device incorporating the 
present invention preferably has an adjustable sensitivity/specificity operating point. 
The sensitivity/specificity operating point is selected by adjusting the threshold in the 
classifier that determines the probability of detection of an acute cardiac ischemic 

15 condition. The sensitivity/specificity operating point may be set by the manufacturer 
of the device, by the purchaser of the device, or by a user of the device. For example, 
if the sensitivity/specificity operating point of the device is set for a high level of 
sensitivity for detection of an acute cardiac ischemic event, the device functions as a 
screening tool that identifies candidates at risk. In a setting where a decision to treat a 

20 patient is being made, the operating point of the device can be readjusted for a high 
level of specificity to confirm the diagnosis of acute cardiac ischemia in the patient. 

Brief Description of the Drawings 
The foregoing aspects and many of the attendant advantages of this invention 
will become -more readily appreciated as the same becomes better understood by 

25 reference to the following detailed description, when taken in conjunction with the 
accompanying drawings, wherein: 

FIGURE 1 illustrates an example of an ECG waveform corresponding to a 

single heartbeat; 

FIGURE 2 depicts a device configured to operate in accordance with the 
30 present invention, with ten electrodes attached to a patient; 

FIGURE 3 is a block diagram of the major components of the device shown in 
FIGURE 2; 

FIGURE 4 is a flow diagram illustrating an acute cardiac ischemia detection 
process conducted by the device shown in FIGURES 2 and 3 which detects acute 
35 cardiac ischemic events by classifying global features derived from a patient's ECG; 
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FIGURE 5 is a flow diagram illustrating a process used to acquire ECG data 
for the acute cardiac ischemia detection process shown in FIGURE 4; 

FIGURE 6 is a flow diagram illustrating a process used to generate 
representative heartbeat data for the acute cardiac ischemia detection process shown 
in FIGURE 4; 

FIGURE 7 is a flow diagram illustrating a process performed during a training 
phase to derive a set of basis vectors used by the acute cardiac ischemia detection 
process shown in FIGURE 4 to detect acute cardiac ischemia; 

FIGURE 8 is a flow diagram illustrating a process performed during the 
training phase to derive a set of classifier parameters used by the acute cardiac 
ischemia detection process shown in FIGURE 4; 

FIGURE 9 is a pictorial diagram of twelve leads of representative heartbeat 
data generated by the acute cardiac ischemia detection process shown in FIGURE 4; 

FIGURE 10 is a pictorial diagram of a concatenated representative heartbeat 
vector x that includes the representative heartbeat data depicted in FIGURE 9 and 
generated by the acute cardiac ischemia detection process shown in FIGURE 4; 

FIGURES 11 A and 1 IB are pictorial diagrams of matrices V and A produced 
during the basis vector derivation process shown in FIGURE 7 by performing a 
Karhunen-Loeve transformation on a correlation matrix R calculated from ECG data 
obtained from a training population of patients; 

FIGURE 12A is a pictorial diagram of twenty projection coefficients ("global 
features") calculated during the acute cardiac ischemia detection process shown in 
FIGURE 4 by projecting a patient's concatenated representative heartbeat vector x, as 
shown in FIGURE 10, onto ten basis vectors spanning acute cardiac ischemia ECG 
subspace and ten basis vectors spanning non-ischemia ECG subspace, the basis 
vectors being determined beforehand as shown in FIGURE 7; 

FIGURE 12B is a pictorial diagram of a concatenated global feature vector f 
that includes the twenty global features depicted in FIGURE 12A; 

FIGURE 13 is a flow diagram illustrating an alternative acute cardiac ischemia 
detection process conducted by the device shown in FIGURES 2 and 3 which detects 
acute cardiac ischemic events by separately classifying global and local features 
derived from a patient and classifying the results of the separate classifications; 

FIGURE 14 is a flow diagram illustrating an alternative acute cardiac ischemia 
detection process conducted by the device shown in FIGURES 2 and 3 which detects 
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acute cardiac ischemic events by classifying a combination of global and local 
features; 

FIGURE 15 is a graph of a typical receiver operating characteristics curve 
depicting sensitivity/specificity tradeoffs which can be implemented by the device 
5 shown in FIGURES 2 and 3 ; and 

FIGURE 16 is a flow diagram illustrating an acute cardiac ischemia detection 
process conducted by the device shown in FIGURES 2 and 3 wherein the 
sensitivity/specificity operating point of the device is adjusted to achieve a desired 
sensitivity/specificity tradeoff as shown in FIGURE 15. 

10 Detailed Description of the Preferred Embodiment 

In regard to the present invention, acute cardiac ischemia includes conditions, 
both chronic and sudden, that result in a deprivation of blood, and hence oxygen, to 
the muscles of the heart, thus requiring urgent treatment if the long-term viability of 
the cardiac muscle cells is to be maintained. While these conditions are reversible at 

15 the time of detection, these conditions will result in permanent damage to the cardiac 
muscle cells if sufficient blood flow is not restored. FIGURE 2 depicts a device 10 
that detects and reports cardiac abnormalities associated with acute cardiac ischemia 
in accordance with the present invention. When the device 10 is attached to a patient 
via a plurality of sensing electrodes (e.g., ten electrodes), the device 10 obtains ECG 

20 data from the patient, automatically evaluates the data, and reports whether the data 
indicates an acute cardiac ischemic condition in the patient. 

The ten electrodes shown in FIGURE 2 are attached to the skin of the 
patient 40 at a variety of locations. First and second sensing electrodes 12 and 14 are 
shown attached to the right and left shoulder areas of the patient 40, respectively. 

25 Third and fourth sensing electrodes 16 and 18 are shown attached to the left and right 
side areas of the patient's torso, respectively, near the patient's legs. A fifth sensing 
electrode 20 is shown attached to the patient's chest area: over the heart. Sixth 
through tenth sensing electrodes 22, 24, 26, 28, and 30 are also spread across the 
chest of the patient 40. 

30 The ten electrodes shown in FIGURE 2 are placed in the standard positions 

for a 12-lead ECG. The first through fourth sensing electrodes 12, 14, 16, and 18 are 
. known as "limb" electrodes. The fifth through tenth sensing electrodes 20, 22, 24, 26, 
28, and 30 are known as "precordial" electrodes. The signals received from the ten 
electrodes are used to produce 12-lead ECG data in a manner well-known to those 

35 skilled in ECG technology. Since the techniques for producing 12-lead ECG data 
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from ten sensmg electrodes is well-known and does not form part of this invention, 
how 12-lead ECG data is produced is not described here. As is conventional, the 
twelve leads are identified as the I, n, HI, aVR, aVL, aVF, VI, V2, V3, V4, V5 and 
V6 leads. ' 

For a user's convenience, the precordial electrodes may be embodied in a 
single patch 42. The patch 42 is designed for quick and easy placement of the 
precordial electrodes on the patient's chest and may be used to reduce the amount of 
cabling needed to connect the precordial electrodes to the device 10. In the event a 
patch 42 is not used, each of the precordial electrodes may be individually attached to 
the patient 40. 

Although FIGURE 2 shows electrode placement in accordance with one 
actual embodiment of the invention, it will be appreciated that the electrodes may be 
placed at alternative locations. For example, the sensing electrodes 12, 14, 16, and 18 
may be placed at the ends of the limbs of the patient 40 as indicated by the open 
circles 32, 34, 36, and 38, respectively. Moreover, it will be appreciated from the 
discussion below that a device formed according to the present invention may use an 
alternative number of electrodes other than ten sensing electrodes (i.e., more or less 
sensing electrodes than shown in FIGURE 2). Regardless of the number of sensing 
electrodes chosen, one or more leads of ECG data are derived from the signals sensed 
by the electrodes. These leads, which may be a subset selected from a larger set of 
leads potentially derived from the ECG signals sensed by the electrodes, are referred 
to herein as the "available leads." 

FIGURE 3 is a block diagram illustrating the major components of the 
device 10 shown in FIGURE 2. ECG signals sensed by the sensing electrodes 
described above (not shown in FIGURE 3) are communicated to a preamplifier 50 via 
lines 52. The preamplifier 50 both amplifies and filters the ECG signals. 
Amplification is required because the strength of the signals sensed by the electrodes 
is generally too low (i.e., in millivolts) to be analyzed by the circuitry of the device. 
The preamplifier 50 may amplify the ECG signals on lines 52 by a factor of 1,000 or 
more. 

The amplified ECG signals are filtered to eliminate noise and other signal 
contaminants. In one actual embodiment of the invention, the filtering includes a 
high-pass filter 54 that attenuates low frequency signals (e.g., frequencies 
below 0.05 Hertz), a low-pass filter 56 that attenuates high frequency signals (e.g., 
frequencies above 150 Hertz), and a notch filter 58 configured to attenuate signals at 
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a particular frequency (e.g., 50 or 60 Hertz, depending on the local line power 
frequency). Alternative embodiments of the present invention may include further 
signal filtering to adapt the device for use in a particular environment. 

The amplified and filtered ECG signals are converted into digital ECG data by 

5 an analog-to-digital (A/D) converter 60. In the embodiment shown in FIGURE 3, the 
ECG signals are multiplexed by the preamplifier 50 and serially communicated to the 
A/D converter 60. Alternatively, the ECG signals may be communicated in parallel 
from the preamplifier 50 to_the A/D converter 60 via separate lines for each of the 
signals (not shown) and multiplexed by the A/D converter. Still further, a separate 

10 A/D converter can be provided for each amplified and filtered ECG signal and the 
output of the multiple A/D converters multiplexed. 

The digital ECG data produced by the A/D converter(s) is applied to a 
processing unit 62 for further processing and evaluation. A memory 64 in 
communication with the processing unit 62 stores the digital ECG data and other data 

15 subsequently generated by the acute cardiac ischemia detection process described in 
more detail below. The memory 64 also stores an acute cardiac ischemia detection 
process 66 in the form of computer program instructions that, when executed by the 
processing unit 62, evaluates the digital ECG data and detects the occurrence of acute 
cardiac ischemia. The memory 64 also stores preprocessed parameters 67 derived 

20 earlier from patients in a training population during a training phase (described 
below). The preprocessed parameters 67 are used by the acute cardiac ischemia 
detection process 66 to evaluate the ECG of a current patient for acute cardiac 
ischemia. 

Although analog signal filters 54, 56, and 58 are shown in FIGURE 3 and 
25 described above, those skilled in the art will appreciate that, alternatively, digital 
filtering after the ECG signals are converted from analog form to digital form can be 
used, if desired. Furthermore, filtering of the ECG data can be performed after the 
data is stored in the memory 64, rather than being performed before storage, as 
shown. 

30 FIGURE 3 also depicts an input device 68 and a display 70 in communication 

with the processing unit 62 for exchanging information with a user of the device. The 
input device 68 allows a user to input information and selectively adjust the operation 
of the device while the display 70 allows the device to report the results of an ECG 
evaluation to the user. The display 70 also permits the device to communicate 

35 instructions to the user. 
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Those of ordinary skill in the art will appreciate that various devices may be 
used to implement the function of the components shown in FIGURES For 
example, the processing unit 62 may be a microprocessor controlled by computer 
program instructions stored in the memory 64. The memory 64 may include 
5 nonvolatile memory in the form of read-only memory (e.g., EPROMs), storage 
memory (e.g., a hard drive), and volatile memory in the form of a random access 
memory (RAM). The input device 68 may include keys, dials, or switches. Similarly 
the display 70 may be a combination of lights or a text display screen, e.g., AMLCD 
LCD, or printer. Audible alerts may also be provided. The construction of suitable 
10 signal amplifiers, filters, and analog-to-digital converters are well-known to those with 
ordinary skill in the art and in many cases are readily available in off-the-shelf devices. 

As part of an analog-to-digital conversion process, ECG signals sensed by the 
electrodes are sampled to obtain discrete voltage values at discrete time intervals 
The rate at which the ECG signals are sampled depends on the configuration of the 
analog-to-digital converter 60. In one actual embodiment of the present invention, 
ECG signals are sampled at a rate of 500 samples per second. Those of ordinary skill 
m the art will appreciate that other sampling rates may be used. 

As will be more readily understood from the discussion below, the device 10 
not only acquires ECG data but also evaluates the data and reports when conditions 
associated with acute cardiac ischemia are detected, including conditions that would 
lead to an acute myocardial infarction or unstable angina. The device 10 extracts and 
evaluates subtle, globaUy-distributed features from a patient's ECG. The extracted 
global features are then classified to determine whether an acute cardiac ischemic 
condition has been detected. In a further embodiment of the present invention, the 
device 10 also derives local features from a patient and classifies both the global' and 
local features to determine whether an acute cardiac ischemic condition has been 
detected. 

FIGURE 4 is a flow diagram illustrating a version of an acute cardiac ischemia 
detection process 66a that determines detects acute cardiac ischemia based on 
• classifying global features. The detection process includes both the extraction of 
global features in blocks 100-114 and the classification of the global features in 
block 116. The detection process begins in a block 100 with the acquisition of ECG 
data from a patient. As illustrated more fully in FIGURE 5, ECG data acquisition 
begins in a block 122 with attaching electrodes, such as sensing electrodes 12, 14, 16, 
18, 20, 22, 24, 26, 28, and 30 to a patient 40. Once the electrodes are attached to a 
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patient, the user of the device 10 initiates ECG data acquisition in a block 124 by 
initiating an input device 68, e.g., by pressing an ANALYZE button on the device. 
Alternatively, the device 10 could automatically initiate ECG data acquisition, e.g., 
upon expiration of a predetermined period of time counted by a timer that is started 
5 upon deployment or activation of the device 10. For a period of time thereafter, the 
voltage signals sensed by the electrodes are amplified in a block 126, filtered in a 
block 128, and converted to digital ECG data in a block 130 as discussed earlier in 
reference to the preamplifier JO and the A/D converter 60 shown in FIGURE 3. The 
ECG data is subsequently stored in a block 132 in the memory 64 of the device 10. 

10 The period of time in which ECG data is acquired in block 100 of FIGURE 4 

is long enough to obtain sufficient, high-quality data representative of one or more 
heartbeats. In one actual embodiment of the present invention, about ten seconds of 
ECG data is acquired. This period of time may be increased or decreased depending 
on factors that affect the quality of the data acquired, including the quality of the 

15 connection between the electrodes and the patient, whether the patient is moving, and 
whether significant electromagnetic noise is present. Accordingly, depending on 
various factors, the period time for ECG acquisition may be increased or decreased 
to, for example, twenty seconds or five seconds, and still be considered an equivalent 
of about ten seconds of data for purposes of the subsequent processing discussed 

20 below. 

After acquiring a patients ECG data in block 100, the device 10 analyzes the 
ECG data and generates representative heartbeat data for each available lead in a 
block 102. Preferably, the representative heartbeat data represents an atrially 
stimulated heartbeat that is common for the lead and is typically characterized by a 

25 high signal-to-noise ratio. If the acquired ECG data is of sufficient high quality, a 
single beat may be used for the representative heartbeat data. In most cases, 
however, it is preferable to combine two or more heartbeats on each lead to generate 
representative heartbeat data for the lead. 

FIGURE 6 illustrates in more detail a process for generating representative 

30 heartbeat data suitable for use in the present invention. In the process shown, 
multiple heartbeats on a lead are used to generate the representative heartbeat data for 
the lead. The device 10 first , locates the heartbeats on each lead in a block 134. In 
that regard, heartbeats may be located by detecting the dominant feature of the QRS 
complexes, i.e., the R waves (see FIGURE 1), of each heartbeat. Because the quality 

35 of the representative heartbeat data is improved by selecting and using only high 
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heartbeat data. For example, heartbeat daa ft at exhibits signiflcam „„. , 

or tue yKS, energy values measured fnr wamni. u 
; i j j • . . f M ' ror exa mple, by a sura of squared values') is 

heartbeat. Of the retained heartbeat data, the middle three heartbeats (middle in a 
temporal sense) are analyzed and two of me three heartbeats that have a jghest 

non-selected third heartbeat .s excluded from further consideration 

1116 r f ned ■»■*— «• to individually compared with the firs, heartbeat 
of foe par selected above. Specifically, a cross-corralation value is determined for 
each of foe QRS combes of me remaining heartbeats with respect to the" 
of the elected pat. Those heartbeat, having a correlation value greater man a 
specmed thread value (e.g., 0.9 ou, of a petfeo, ,.0) are retained. 
havtng a lower correlation value are excluded in a block 140 from farther 
onto, for generating a represent heartbeat The heartbeats reamed are 
used in generating a representative heartbeat for the lead being analyzed 

To generate a representative heartbeat for the lead, the waveforms of foe 
retained heartbeats are aligned in a block 142 by matching up wavefortn featutes such 
* fo R wave peak, preferably wording to me maximum oross- M „elation of data 
samp as m foe reamed heartbeaa. Further, a baseline offset oalculated for each ft 
. subtracad m a block 144 from foe daa samp.es in eaoh reamed heartbeat £ 
basckne offset ts foe amouu, foa, foe daa samples in a heartbeat exceed a selected 
efeence level. In au actual embodiment of foe invention, foe basehne offset is 
estimated as foe mean of 48 milliseconds of ECG data taken from an 88 millisecond 
mark to a40 mtihsecond mark prior to foe QRS peak sample (e.g., if foe QRS peak 
™ . located a, foe 100 millisecond mark, foe estimated baseline offse, was foe mean 
of the samples from foe 12 millisecond mark to foe 60 millisecond mark) I, is 
apprecaad that more or leas ECG daa may be used, or foe timing of foe ECG 
segment used may be shifted, in calculating a baseline offset 

Next, in a block 146, a set of values repreaenative of an average heartbeat for 
each lead is computed from foe retained heartbeat data. In an actual embodiment of 
foe invention, foe repreaenative heartbeat data is generated by calculating a sample- 
by-sample average of the retained heartbeats. Nevertheless, it is understood foa, a 
wide vanety of mafoematical operations may be used instead of averaging For 
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example, representative heartbeat data may be calculated from a sample-by-sampie 
mean, mode, weighted mean (using weighing coefficients), or trimmed mean (i.e., 
calculate a mean from a set of values that excludes values that are too large or too 
small), or median. 

5 Before the resulting representative heartbeat data is subjected to further 

evaluation, it may be advantageous to reduce the amount of data being evaluated. In 
that regard, data samples in a segment of the representative heartbeiat are selected in a 
block 148 to represent the.ECG information of particular interest. In one actual 
embodiment of the present invention, a 480 millisecond segment of representative 

10 heartbeat data for each available lead is selected beginning at an 8 millisecond mark 
before each representative beat's QRS peak. Data outside the selected segment is 
excluded from the further processing and evaluation performed. In reducing the 
amount of data being evaluated, it is appreciated that a segment of ECG data longer 
or shorter than 480 milliseconds may be selected. The beginning mark may also be 

1 5 adjusted earlier or later in the sequence of ECG data. 

While the process shown in FIGURE 6 is used to generate representative 
heartbeat data in one actual embodiment of the invention, it should be understood that 
other ones of the many known procedures for generating representative heartbeats 
can be used. For example, the methods of representative heartbeat generation used in 

20 the LEFEPAK® 11 or LIFEPAK® 12 defibrillators manufactured by Physio-Control 
Corporation of Redmond, Washington, may be used, if desired. 

Returning again to FIGURE 4, once representative heartbeat data is generated 
for each of the available leads in block 102, in accordance with the present invention, 
the data is subjected to a global classification evaluation in blocks 104-116. As will 

25 be better understood from the discussion below, a global classification evaluation 
involves concatenating the representative heartbeat data produced on the available 
leads, extracting global features from the concatenated representative heartbeat data, 
and determining whether the global features are in a class indicative of acute cardiac 
ischemia. More specifically, as will be discussed more fully below, a global 

30 classification evaluation involves mathematically projecting a concatenated vector of 
representative heartbeat data onto predetermined basis vectors that define an acute, 
cardiac ischemic ECG subspace and onto predetermined basis vectors that define a 
non-ischemic ECG subspace. The resulting projection coefficients are the global 
features that are classified to determine whether an acute cardiac ischemic condition is 
35 detected. 
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an actual embodiment of the invention k block I48 <Fra u»B 6), i, is again noted to 

any segment lengti, (preferably containing, a, leaat, the QRST portion! of " 
may be used. Nevertheless, ae ™h be better understood front the discussion bel!w 
fixed number of samples in a represent heartbeat mua be set aprior, during' 

Ptee for fc to Permed global damnation evaluation of 
the present tnvention. The length of tin. basis vectors detetmineti during the framing 
phase and later used in the-global classification evaluation depends on both the 
number of data samples in each representative heartbeat and the number of available 
eads use from the ECGs of a training population of patients in the preprocessing 
training phase. F s 

Prior to discussing the global classification evaluation illustrated in 
blocks 104-116 of FIGURE 4, a brief discussion of son* concepts underlying the 
evaluation is provided, A description of the preprocessing training phase required to 
15 derive the basis vectors is also provided. 

The global classification evaluation shown in FIGURE 4 is based on the 
concept that a series of numbers can be viewed as coefficients of a vector that defines 
a point in a multidimensional signal space. In the case of the present invention ' this 
means that a series of ECG data points in a patient's representative heartbeat may be 
0 viewed as coefficients of a vector defining a point in a multidimensional ECG signal 
space. The number of dimensions in the ECG signal space is determined by the 
number of data points in the patient's representative heartbeat data. In a similar 
manner, a series of global features derived from a patient's ECG according to the 
present invention may be viewed as coefficients of a vector defining a point in a 
5 multidimensional feature space. The number of dimensions in the feature space is 
determined by the number of global features derived from the patient 1 s ECG It 
should be understood that the ECG signal space and the feature space are completely 
different conceptual spaces, though in regard to this invention they are related by 
virtue of the projecting operation described later in reference to blocks 1 1 0 and 1 12 
» of FIGURE 4. 

For purposes of discussion herein, patients in a training population of patients 
known to have acute cardiac ischemia are referred to as ischemic training patients or 
an ischemic training population. Likewise, patients in a training population of patients 
known to not have acute cardiac ischemia are referred to as non-ischemic training 
patients or a non-ischemic training population. Sets of global features derived from 
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the ECG data of ischemic training patients may be plotted as points defining an 
"ischemic" region of a feature space. In a similar fashion, sets of global features 
derived from the ECG data of non-ischemic training patients may be plotted as points 
defining a "non-ischemic" region of the feature space. According to the present 
5 invention, if a set of global features derived from the ECG data of a patient under 
current evaluation defines a point closer to the ischemic region than the non-ischemic 
region, a diagnosis of acute cardiac ischemia is reported. The threshold of "closeness" 
for making such a diagnosis, is adjustable, thus providing the device with an adjustable 
sensitivity/specificity tradeoff, as will be discussed in more detail below. 

10 Within the ECG signal space (which, as noted, is different than the feature 

space), ECG data obtained from ischemic training patients is viewed as defining an 
ischemic ECG subspace. Likewise, ECG data obtained from non-ischemic training 
patients is viewed as defining a non-ischemic ECG subspace. The ischemic and 
non-ischemic ECG subspaces may be succinctly characterized by mathematical basis 

15 vectors derived in a training phase (discussed below). The basis vectors are stored as 
preprocessed parameters 67 in memory 64 of the device 10. As will be discussed later 
in greater detail, the basis vectors are used in blocks 1 10 and 112 of FIGURE 4 to 
derive the projection coefficients (i.e., global features) of the patient under current 
evaluation. 

20 Before deriving a set of basis vectors for characterizing each of the ischemic 

and non-ischemic ECG subspaces, the quality of this characterization may be 
enhanced by dividing the ischemic and non-ischemic training populations into smaller 
groups according to a selected characteristic. Sets of basis vectors (ischemic and 
non-ischemic) are then derived for each group from the ECG data of the training 

25 patients in the group. However, dividing the training populations into groups is not 
required. 

In one embodiment of the invention, the training populations are divided into 
groups according to locations of ischemic conditions (e.g., anterior, inferior, and 
other). The selected characteristic used to divide the training populations is the 
30 identity of the lead with the greatest ST elevation. It has been found that the identity 
of the lead with the greatest ST elevation is generally indicative of the location of an 
acute cardiac ischemic condition. Accordingly, each patient in the ischemic and 
non-ischemic training populations is assigned to an "inferior/' "anterior," or "other" 
group based on the patient's ST elevation measures. 
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These leads are depicted generally in FIGURE 9 in which dots are used to indicate a 
series of numbers that, in this regard, forms the representative heartbeat data. 

Assuming for the sake of discussion that the training populations are divided 
into groups according to locations of acute cardiac ischemia as described above, the 
5 ECG of each patient in the training populations is categorized into a group in a 
block 154 according to the ST elevation measures calculated for each patient. Next, 
in a block 156 the representative heartbeat data generated for each patient in each 
group is concatenated to form a representative heartbeat vector "x". A concatenated 
representative heartbeat vector x includes the representative heartbeat data of a first 

10 lead (e.g., lead I), immediately followed by the representative heartbeat data of a 
second lead (e.g., lead II), a third lead (e.g., lead m), a fourth lead (e.g., lead aVR), 
and so on. A concatenated representative heartbeat vector x is depicted generally in 
FIGURE 10, with vertical dashed partitions between the lead data shown for 
illustrative purposes only. It is to be understood that alternative embodiments of the 

15 invention may use a different (and/or smaller) combination of leads than that shown 
inFIGURES9andlO. 

In an embodiment of the invention where twelve leads of ECG data are 
acquired at a sampling rate of 500 samples per second, and a 480 millisecond segment 
of data per representative heartbeat is used, each lead (e.g., the leads shown in 

20 FIGURE 9) includes 240 samples of heartbeat data, thus producing a concatenated 
representative heartbeat vector x (as shown in FIGURE 10) having 2,880 samples of 
heartbeat data. A concatenated representative heartbeat vector x is thus obtained for 
each of the patients in the respective training populations. 

For the ischemic and non-ischemic sets of patients in each group, the patients' 

25 concatenated heartbeat vectors x are combined together in a block 158 in FIGURE 7 
to calculate correlation matrices R^ and R^ using the following general equation: 

i = l 

The values in matrix R of Equation (1) are normalized according to the total 
number of patients "N" whose concatenated heartbeat vectors x were used. Thus, 
30 using Equation (1), a correlation matrix R^c is calculated for the ischemic patients in 
each group, and a correlation matrix R non is calculated for the non-ischemic patients 
in each group. 
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vectors and a set of non-ischemic basis vectors) may be used to represent the ECG 
subspaces of the ischemic population and the non-ischemic population as a whole. 
The selected vectors are then stored in a block 163 as preprocessed parameters 67 in 
memory 64 for later use in deriving global features from the ECG of a patient 
5 undergoing evaluation. 

Returning now to FIGURE 4, given a defined number of groups and a paired 
set of ischemic and non-ischemic basis vectors associated with each group, the ECG 
of a patient under current evaluation is categorized in a block 104 into the group to 
which it best pertains. As part of categorizing the patient's ECG, composite 

10 ST elevation measures uj, U2, and U3 are calculated for the current patient as they are 
described above for the training populations. The current patient ! s ECG is 
categorized into the "inferior" group if the measure uj is greatest. If the measure u 3 is 
greatest, the patient's ECG is categorized into the "anterior" group. If the measure u 2 
is greatest, the patient's ECG is categorized into the "other" group. 

15 Once a patient's ECG is categorized into a particular group in block 104, the 

paired set of basis vectors associated with the group are selected in a block 106 from 
the preprocessed parameters 67 in memory 64 for use in evaluating the patient's ECG. 
The patient's representative heartbeat data for each of the available leads is then 
concatenated in a block 108 to form a concatenated heartbeat vector "x", in a manner 

20 as described earlier for the patients in the training populations. Next, in blocks 1 10 
and 1 12 of FIGURE 4, the patient's concatenated heartbeat vector x is mathematically 
projected onto the basis vectors selected in block 106. Broadly stated, the projecting 
operation (described in more detail below) results in a number of projection 
coefficients that are used as global features of the patient's ECG. The number of 

25 global features extracted from the patient's ECG corresponds to the number of basis 
vectors used in the projecting operation. 

More specifically, the projecting operation involves projecting the patient's 
concatenated heartbeat vector x onto the basis vectors defining an acute cardiac 
ischemic ECG subspace in block 110 by computing an inner product of the vector x 

30 with each of the ischemic basis vectors. The patient's concatenated heartbeat vector x 
is also projected onto the basis vectors defining the corresponding non-ischemic ECG 
subspace in block 1 12 by computing an inner product of the vector x with each of the 
non-ischemic basis vectors. If, for example, ten basis vectors are used to characterize 
each of the ischemic and non-ischemic ECG subspaces, the projecting operation 

35 results in a total of twenty scalar projection coefficients that are used as global 
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training patient are concatenated to form a global feature vector f in the manner as 
earlier described and generally depicted in FIGURE 12B. 

Mean feature vectors and m non are then calculated in a block 167 for 
each group by computing a feature-by-feature mean of the concatenated global 
5 feature vectors f produced from the respective ischemic and non-ischemic patients in 
each group. In a configuration where ten basis vectors are used to characterize each 
of the ischemic and non-ischemic ECG subspaces, the concatenated global feature 
vector f for the patients in_the training populations, and thus the mean feature 
vectors m^ and m non for each group, will include twenty values. 
10 In addition to calculating mean feature vectors m^ and m non in block 167, 

co variance matrices Cj sc and C non are calculated in block 168 for the respective 
ischemic and non-ischemic training patients in each group. Generally speaking, a 
covariance matrix C contains the statistical covariance of the global feature vectors f 
of a set of N number of patients, and is calculated as follows: 

15 C = l Zfcfi-m^i-ffi) 7 ) (3) 

The covariance matrix C provides an indication of how widely dispersed the 
global features of the set of patients are from the calculated mean global features for 
the same set of patients. The mean feature vector m and the covariance matrix C thus 
calculated for each of the ischemic and non-ischemic sets of patients in each group are 

20 stored in a block 169 as preprocessed parameters 67 in memory 64 (FIGURE 3) for 
later use in classifying the current patient's global features. 

Returning to FIGURE 4, according to one implementation of the invention, 
the Gaussian classifier in block 116 evaluates the current patient's global feature 
vector f with respect to the mean global feature vectors m^ and m non , and covariance 

25 matrices Cj sc and according to the following equation: 

d = ^/(f~m) T C- 1 (f-m) (4) 

The quantity ,f d H computed from Equation (4) reflects the distance between 
the patient's global feature vector f and a mean global feature vector m weighted by a 
covariance matrix C. Accordingly, a quantity d^ reflects the distance between a 
30 patient ! s global feature vector f and the mean global feature vector m^ weighted by a 
covariance matrix for the acute cardiac ischemia ECG subspace. Similarly, a 
quantity d non reflects the distance between a patient's global feature vector f and the 
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if the training populations are divided into groups according to ischemia location, the 
reported outcome may also identify the location of the ischemic condition (if acute 
cardiac ischemia is detected) based on the group into which the patient's ECG was 
categorized. Knowing whether a detected ischemic condition is at an inferior, 
5 anterior, or other location may assist a caregiver in treating the ischemic condition. 

Although the acute cardiac ischemia detection process 66a described above 
involves projecting a patient's concatenated heartbeat vector x onto basis vectors that 
collectively define both an acute cardiac ischemic ECG subspace and a non-ischemic 
ECG subspace, it should be understood that, alternatively, the patient's concatenated 
10 heartbeat vector x may be projected onto one or more basis vectors that define only 
an acute cardiac ischemic ECG subspace, i.e., without projection onto any basis 
vectors that define a non-ischemic ECG subspace. In that regard, only "ischemic" 
projection coefficients (i.e., ischemic global features) are produced and classified. 

A classifier for classifying only ischemic global features does not need to be 
15 any different in structure than a classifier that classifies both ischemic and non- 
ischemic global features, as described above in reference to block 116 (FIGURE 4). 
The only difference is the number of global features used in the classification and the 
training of the classifier performed beforehand in a training phase. 

In the training phase, the basis vectors defining an acute cardiac ischemic ECG 
20 subspace are derived from the ischemic training population as described earlier in 
reference to FIGURE 7. Then, for each training patient in both the ischemic and non- 
ischemic training populations, ischemic global features are calculated using the 
derived ischemic basis vectors in a manner as described earlier in reference to 
block 164 in FIGURE 8. The ischemic global features for each patient (ischemic and 
25 non-ischemic) are concatenated into a global feature vector f as described earlier in 
reference to block 166. 

Mean feature vectors m^ and m non and covariance matrices Qsc and Cn 0n are 
then calculated in a manner as described in reference to blocks 167 and 168. The 
vectors m^ and m non , and matrices Cue and C n0 u, are stored as preprocessed 
30 parameters 67 in the memory 64 for later use in classifying the ischemic global 
features of the current patient as discussed in reference to block 116. 

Alternatively, the patient's concatenated representative heartbeat vector x may 
be projected onto one or more basis vectors that define only a non-ischemic ECG 
subspace, i.e., without projection onto any basis vectors that define an acute cardiac 
35 ischemic ECG subspace. In that regard, only "non-ischemic" projection coefficients 
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To determine the set of local features, one or more local morphological 
measures are calculated from the representative heartbeats in block 184. Patient 
clinical information, such as the age and sex of the patient, may also be obtained in 
block 186 from the user of the device 10. With regard to sex information, a "1" may 

5 be used for males and a "0" may be used for females. If clinical information is not 
entered or available at this time, default values may be used. While clinical 
information may improve accuracy in the detection of acute cardiac ischemia, the 
ECG classification of the pre§eiit invention can be performed without this information. 
The values of the one or more local morphological measures and clinical information 

10 are concatenated to form a local feature vector that is input into a local feature 
classifier in a block 188 that computes a local classification, statistic. The local 
classification statistic is not compared with a decision threshold at this time, but 
instead is also provided to the combiner in block 190. 

The local feature classifier in block 188 may be either a statistical classifier or 

15 a heuristic classifier. A suitable statistical classifier for this aspect of the invention 
may be one similar in form to the Gaussian classifier used for classifying global 
features and described above. In that regard, the local feature classifier in block 188 
is trained in a similar fashion using selected local morphological measures and clinical 
information as features : derived from : both ischemic and non-ischemic training 

20 populations. The selected local features for each training population, or group within 
a population, are concatenated to form local feature vectors from which a mean local 
feature vector is calculated. The corresponding local features of the patient are 
evaluated, using Equation (4), against the mean local feature vectors. The result of 
this -evaluation is a local classification statistic that is provided to the combiner in 

25 block 1 90 as noted above. 

Alternatively, a heuristic classifier may be used. A heuristic classifier tries to 
mimic procedures used by an expert (i.e., a cardiologist) to evaluate an ECG and 
report a diagnosis. A heuristic classifier compares the local features to expert- 
determined thresholds. For an expanded description of a suitable heuristic classifier, 

30 see G. Wagner, Marriott's Practical Electrophysiology, 9th Ed. (1994), published by 
Williams & Wilkins, Baltimore, which is incorporated herein by reference. The result 
of the evaluation is a local classification statistic that is supplied to the combiner in 
block 190. 

As another alternative, or in addition to directly, inputting a patient's local 
35 morphological measures and clinical information as local features into a local feature 



WO 99/55226 
« 



WO 99/55226 



-26- 



PCT/US99/09311 



10 



15 



20 



25 



_ ^ CIinical inforaati0fl may be caIcdated . ^ 

dasher. In one actual embodiment of the invention two different nrln 
used for ^^ 

in rri r duce a pre,iminaiy predicti ° Q ° f ^ 

schemia is pre sat The Iogistxc regression equation is derived during a preprocessing 
^ngphase(,g., the tranmgphase described earlier in which the IJsZmZ 
c orner parameters, i.e., the processed parameters 67, are derived) accortg t7 
logistic regression model defined by the following equation. 



30 



K 

a,+ 2 a l x i ( 5 ) 



In Equatanffl, V denote the ^ability of daeoted 
■schemta, a. . a coated consunr. , denotes the ce.cnia.ed id. j£l 

ItT o ( ' e '" h , T, Pr0blbility ° f h 1 -hemic 

ST? J 7 r ^ Pr ° baM,y * ^ »— °^ 

mm m the respect™ taming portions, .he regression coefficients 0 f ft. 

(TIGURE3) ftr later use m producfag , rf ^ 

tschetmamthepadenttindercurrentevalnatioa 

otiwl P ° PUla,i0 ° S *" iM0 «™f «*«**. M»ior and 

«ta) as deacnbed earlier, a logistic regression equation is derived for each oft 

P*. b regard, ECG feantres and leads fcnowu to be associated v^al^ 

.schennc events (e. g ., ft e lMds used „ 
STdevatan measure n 3 discnssed above) are se.ec.ed for use in derivfa! ft 

nfenor or ofter location (e.g„ th, ,eada „sed fa calcufating fte composite 
ST e, m ,o„ meaanres n, or u 2 discussed above) are selected for u* fa derivta. an 
mfenor- and an -ofter- , ogistic regrM , 0 „ ^ ^ = - 

denved for eacb of fbe groups „ ao r=d as preprocessed paramours 67 fa 2 
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memory 64 for later use in calculating a probability of detected acute cardiac ischemia 
in the patient under current evaluation. 

For the patient under current evaluation, the patient's local features are input 
into the logistic regression equation derived during the training phase described 

5 above. More specifically, the patient's local features are weighted by the derived 
regression coefficients and combined as shown in Equation (5) to produce an output 
"p." The output "p" (i.e., the probability of detected acute cardiac ischemia) is used 
as a composite local feature to be classified in the local feature classifier in block 188. 
Where a logistic regression equation is derived for each of an anterior, inferior, and 

10 other group, a probability of detected acute cardiac ischemia in the current patient 
(i.e., a composite local feature) is calculated with respect to each of the groups. 
Thus, in that regard, an "anterior" composite local feature, and "inferior" composite 
local feature, and an "other" composite local feature are calculated for the current 
patient and input into the local feature classifier. 

15 Alternatively, rather than direcdy inputting the composite local features into 

the local feature classifier in block 188, it may be advantageous to first dichotomize 
the composite local features. For instance, if the composite local feature sensitive to 
anterior acute cardiac ischemia exceeds a preselected threshold, the composite local 
feature is assigned a value :of "1" (indicating predicted anterior acute cardiac 

20 ischemia). If the composite local feature does not exceed the threshold, it is assigned 
a. value of "0" (indicating predicted non-ischemia). Similarly, dichotomized outputs 
"1" and "0" may be determined for the composite local features sensitive to inferior 
and other locations of acute cardiac ischemia. The thresholds used to dichotomize the 
composite local features are selected based on prediction patterns observed in the 

25 training populations during the training phase. The dichotomized composite local 
features are then input into the local feature classifier in block 188. For additional 
description of constructing and implementing logistic regression models, see 
D. Hosmer and S. Lemeshow, Applied Logistic Regression (1989), John Wiley & 
Sons, New York, incorporated by reference herein. 

30 Another procedure for creating a composite local feature involves calculating 

a Mahalanobis distance. A Mahalanobis distance "d" is calculated according to 
Equation (4) and measures the distance between the patient's set of local features and 
a representative . set of local features derived from a training population. A 
Mahalanobis distance <!« and/or d non , calculated using representative local features 
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In implementations of the invention where Err. ' . 

(e.g., anterior, inferior, and other) a MahZT „ Ca ' eg ° rizCd **> 
5 patienfs iocai^eawth^o ^ofl " " ' 

Using a neasest-neighbo, apprTh 1^. * * "* n " 1 - «*- groups, 

a group (either antlr ^0^ ^ ^ " *" *■«■* -» 

dichotomized oontpoI o^Z TtJ" - i,ldUdi " 8 ° M « — 

15 Equation (5) to detetmine acute cardiac IT. * ™ m rrfffsnce ,0 

composite local features mTnLT , T ° m « «he 

a ^stic region ^^Xle^ 0 "^ * 
composite local features derived fern patients^ ,!? ' ? ^ , " e 
known outcome of iachemia or «*£Th £ ^ «*> 
20 coefficients are stored as oreoroJs^ ^^ 0115 ' The 

for later use in dJ^T* f » * —T c*4 (FIGURE 3) 

under current cvahtation The p ^ "** "*■* h 

the local feaoure corner in t^!, m" ^ " 

•o the combiner in bio* J 90 for l!t T " sudstic ** " supplied 

25 Caseation statistic. " C ° mbim,iOT «» «* *M 

The combiner in block 190 k * ^ 

Katies torn both me g,obal attd « * 

188 resDectiveM n,. . u- • classifiers (i.e., from blocks 182 and 

statistic^' £ TTZT^ 1 *« - * 

» classifier The coZ- , ? * dWi&r 0r ^ ^^ate 

produceaco^rr^r:!?*" ,o i - 

* -mbined dassmcation ' 
classification statistics are classified „ M ' 0Cal aad « loW 

acute cardiac ischemia ^1^7^ ~» 
5 less than the threshold Y thlt 7,,, 1 classification statistic is 

threshold t , the local and global classification statistics are classified as 
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belonging to a non-ischemic class of patients. Suitable classifiers for use in the 
combiner in block 190 are discussed by R. Duda and P. Hart in Pattern Classification 
and Scene Analysis, referenced above. The result of the classification made by the 
combiner (i.e., whether or not acute cardiac ischemia is present) is then reported to 
5 the user in a block 192. 

FIGURE 14 is a flow diagram that illustrates another alternative version of an 
acute cardiac ischemia detection process 66c formed in accordance with the present 
invention in which both global and local features are evaluated. In FIGURE 14, a 
single classifier is used in a block 196 in place of the separate global and local 

10 classifiers (blocks 182 and 188, respectively) and the combiner (block 190) described 
in reference to FIGURE 13. The single classifier in block 196 receives both the global 
features computed in block 180 and the local features calculated in block 184 as 
features to be jointly classified. Patient clinical information may also be input into the 
classifier in block 196 as additional local features (though not shown in FIGURE 14). 

15 Alternatively, or in addition to the local features discussed above, composite local 
features may be calculated and provided to the classifier in block 196 along with the 
global features.- 

The classifier in block 196 may be a statistical classifier having a form similar 
to the classifiers described, above. For example, the classifier may be a Gaussian 

20 classifier previously trained during a preprocessing training phase using corresponding 
sets of global and local features derived from patients in the ischemic and non- 
ischemic training populations. The current patient's global features and local features 
are concatenated into a single global/local feature vector. The classifier then 
evaluates the patient's combined global/local feature vector with respect to calculated 

25 representative global/local feature vectors derived from the ischemic and non-ischemic 
populations to produce a classification global/local statistic. The global/local 
classification statistic is then compared with a selected threshold to produce a 
diagnosis of whether an acute cardiac ischemic condition is present. In an actual 
embodiment of the present invention, a logistic regression classifier is used to evaluate 

30 the combined global/local feature set. For additional detail on a logistic regression 
classifier, see D. Hosmer and S. Lemeshow, Applied Logistic Regression (1989), 
referenced above. Depending on the outcome of the evaluation performed by the 
classifier in block 196, a report indicating detection of acute cardiac ischemia is 
produced in block 1 92. 
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It should be understood that a classifier Utvm^ii • , 
software routine. In reference to FIGURE^? * ^ ^plemented as a computer 

sensitivity/specificity tradeoff i 5 fflJLd by "er t^T T" ' 
(ROC) cu™ 200 as shown in FIGURE ,5. Cl io 
sensitivity, or fraction of true positives detected JL 

nf ■ i » . pMraw aetect ed, and the x-ajas represents the quantitv 

of 1- spectfictty", or fraction of to positives detected. Values on both ITar! 
expressed as percentages. Thus, as shown in FIGURE 15 . a., . , 

to be more speofic is typically tea se^^. If ^ ' * „ , 

In tit. present utvenUon, the point « wtich a detection device , 0 operateswtftsROC 
curve ntay be adjusted by varying the decision threshofd V IZZ^Z 
caIou 1 a te fdas sifi cations,ati S tic(d e scribedearher)i S co m p^ ' 
«. V "7" f a,0dimeM ° f i""-^ the decision threshold, and hence the 

r . * .tetnative,, the threshold LThe ^^7^^ 
dm. accordtng to the purchase* neetis by adjusting approprie variable It 
asfiware or by setttng an internal dial or switch that is reati by the softwtne ThI 
dcv.ce ,0 »y nsrtiter be contigured to receive a user input (4 an ^ta, m 
«*k or key input, tha, ^actively adjust tite threshold u J by tite £Z a^ 
thus adjusts the sensitivity/speciacity operating point of tite device 

according 10 " "-WW*** tradeoff 

to tit. resent tnvennon May be uaed for screening patients whh possible 
«| c^dtac tschenna as wall aa contoing a detected ischemic condition for 
— purposes. FIGURE ,6 is a flow diagram iterating . acut. cardiac 
- he™, detectronproccss 220 for screening and confining adetected ^ ^ 
.*« condtraon In a bloch 222, a desired sensitivhy/speciftcity operating point for 
tit. devtc. » se^ed etther by the n.nufacu.rer. purchaser, or user, fa tita. re^rd, to 
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screen patients for possible acute cardiac ischemia, it is generally preferred that the 
sensitivity/specificity operating point of the device be set for a high level of sensitivity. 

In a block 224, the device 10 acquires ECG data and, in a block 226, the 
device evaluates the ECG data to detect acute cardiac ischemia as described by the 

5 acute cardiac ischemia detection processes 66a, 66b, or 66c shown in FIGURES 4, 
13, or 14. In a block 228, the device 10 reports the results of the acute cardiac 
ischemia detection process via a user display. In a block 230, if acute cardiac 
ischemia is detected, the user of the device is given opportunity to adjust the 
sensitivity/specificity operating point to a higher level of specificity for a second 

1 0 evaluation of the patient. 

After adjusting the sensitivity/specificity operating point in block 230, the 
device 10 performs a second evaluation of the patient by returning to block 224 to 
acquire ECG data from the patient: The patient's ECG data is analyzed and evaluated 
in block 226 using the adjusted sensitivity/specificity operating point to confirm 

15 whether an acute cardiac ischemic condition has been detected. The result of the 
second evaluation is reported in block 228. It is appreciated that additional 
evaluations of the patient may be performed as desired. 

Although FIGURE 16 illustrates a second evaluation being performed on 
newly acquired ECG data (i.e., by returning to block 224 from block 230), it should 

20 be understood that a second and additional evaluations may be performed on the 
original set of ECG data acquired from the patient in the first evaluation. In that 
regard, rather than returning to block 224 from block 230, the device 10 would return 
to block 226 to reevaluate the original ECG data using the adjusted 
sensitivity/specificity operating point. 

25 While various embodiments of the invention have been illustrated and 

described herein, it is appreciated that changes may be made without departing from 
the spirit and scope of the invention. For example, rather than concatenating the 
representative heartbeat data on each of the leads and projecting the concatenated 
heartbeat vector onto the basis vectors, each of the particular leads may be 

30 individually projected onto a set of basis vectors derived for the particular lead and 
the resulting projection coefficients may be used as global features. It should also be 
understood that with sufficient human expert evaluation of global features, classifying 
the set of global features may be performed heuristically. Furthermore, when 
reporting whether an acute cardiac ischemic condition is detected, a device 

35 constructed according to the invention may report an outcome in a range of outcomes 
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the range of values ,s reported to the nser. It is intended, therefore, that the », of 
the tnvenuon be deternnned from the dahna to Mow and enuivale„ ts 
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The embodiments of the invention in which an exclusive property or privilege 
is claimed are defined as follows: 

1 . A method of detecting and reporting a condition associated with acute 
cardiac ischemia in a patient, comprising: 

(a) obtaining one or more leads of ECG data from a patient; 

(b) forming a vector of heartbeat data derived from the patients 
ECG data; - 

(c) producing a set of global features by projecting the vector of 
heartbeat data onto one or more basis vectors that define an acute cardiac ischemic 
ECG subspace or a non-ischemic ECG subspace; 

(d) classifying the set of global features to determine whether the 
global features are indicative of an acute cardiac ischemic condition; and 

(e) reporting whether an acute cardiac ischemic condition is 

detected. 

2. The method of Claim 1, wherein forming the vector of heartbeat data 
includes: 

(a) analyzing the one or . more leads of ECG data to identify one or 
more heartbeats; 

(b) generating representative heartbeat data for each lead; and 

(c) concatenating the representative heartbeat data for each lead to 
form the vector of heartbeat data. 

3. The method of Claim 1, wherein producing the set of global features 
includes: 

(a) calculating an inner product of the vector of heartbeat data and 
one or more basis vectors that define an acute cardiac ischemic ECG subspace to 
produce a corresponding number of ischemic condition projection coefficients; 

(b) calculating an inner product of the vector of heartbeat data and 
one or more basis vectors that define a non-ischemic ECG subspace to produce a 
corresponding number of non-ischemic condition projection coefficients; and 

(c) using the ischemic condition projection coefficients and the 
non-ischemic condition projection coefficients as the set of global features. 



4. The method of Claim 1 , further comprising: 
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(a) defining a plurality of groups wherein each group has basis 
vectors associated therewith that define an acute cardiac ischemic ECG subspa ce Id 
a non-ischemic ECG subspace; sucspace and 

of groups; and ® * ^ *» ht ° ' -* « «* 

„„. . , . W . UShs ,he basis v « ors ° f «"» group tato which the patient's 

S£r ci,es mhe bMis veaors ° n '° wbkh "» ™ - 

(a) selecting a local feature derived from the patient- and 
sheeted iocal fie. ^ ^ ^ ** to ' ^ b »< - - 

6. The method of Claim 4, further comprising: 

Unrrf f ^ defi ^ geaChgrOUP ° fthe P ,Uralit y of Wtocorrespondto 
a location of an acute cardiac ischemic condition; and - 

(b) if an acute cardiac ischemic condition is detected, reporting the 
locauon of the acute cardkc ^ ^ P^i£ 

which the patient's ECG data is categorized. 

7. The method of Claim 6, wherein the location of an acute cardiac 
aschemic condition is determined based on an ST elevation measured on a lead and 
wherein categorizing the patient's ECG data into a group includes: ' 

nh. 9 - „ * ? • ^ ST deVati0Q ° f ° ne or raore of the leads 

obtained from the patient; 

calculated- ^ f0 ™ S ° f ^ **** ** ^ ST elevation was 

(c) calculating a composite ST elevation for each subgroup by 
calculating a mathematical combination of the ST elevation of the leads in each 
subgroup; and 

(<0 categorizing the patient's ECG into a group according to the 
subgroup whose composite ST elevation is greatest. 

8. The method of Claim 1, further comprising: 
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(a) using a Karhunen-Loeve transformation to calculate a first set 
of basis vectors that define an acute cardiac ischemic ECG subspace and a second set 
of basis vectors that define a non-ischemic ECG subspace; 

(b) selecting one or more basis vectors from the first set of basis 
vectors and one or more basis vectors from the second set of basis vectors as the basis 
vectors onto which the vector of heartbeat data is projected. 

9. The method of Claim 1, wherein classifying the set of global features 
includes: 

(a) concatenating the set of global features to form a global feature 

vector; 

(b) producing a classification statistic by evaluating the global 
feature vector relative to a representative global feature vector derived from a training 
population; and 

(c) comparing the classification statistic with a threshold to 
determine whether an acute cardiac ischemic condition is detected. 

10. The method of Claim 9, wherein a Gaussian classifier is used to 
evaluate the global- feature . vector relative to the representative global feature vector 
and produce the classification statistic. 

11. The method of Claim 9, wherein the threshold is selected in 
accordance with a desired sensitivity/specificity tradeoff. 

12. The method of Claim 1, wherein if the one or more basic vectors 
define only an acute cardiac ischemic ECG subspace, or if the one or more basis 
vectors define only a non-ischemic ECG subspace, then classifying the set of global 
features includes: . 

(a) concatenating the set of global features to form a global feature 

vector; 

(b) producing a classification statistic based on a calculated 
distance between the global feature vector and a representative global feature vector 
derived from a training population; and 

(c) comparing the classification statistic with a threshold to 
determine whether an acute cardiac ischemic condition is detected. 
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13 A device for detecting ^ reporti a condition 
cardiac ischemia in a patient, comprising: 

(a) aplurality of electrodes for sensing patient ECG signals- 

(b) a processing unit in communication with the plurality of 
e ectrodes for evaluating heartbeat data derived from the sensed patient ECG ^al 

rXtT 7 Pr0dUCin§ 3 " ° f gl ° bal ^ by **** * -tor oft 

ECG 1 ° r , m0re ^ V6Ct0rS ** **" - acute ischemic 

ECG subspace or a non-ischemic ECG subspace; 

. f i u , , 3 ClaSSifier CB " CUtod by ^ processin S unit ^ evaluating the 

set of global features and determining whether the global features are indicative 0 f 1 
acute cardiac ischemic condition; and 

r . rt - u " ^ ° UtpUt ^ coramunicati °n with the processing unit for 

reporting whether an acute cardiac ischemic condition is detected. 

14. The device of Claim 13, wherein the classifier executed by the 

(a) a second classifier executed by the processing unT for 
classing a set of local features derived from the patient to produce a second output; 

a, <i « • (b l a f rd dassifier " wcuted by the processin s UQit for r «*™S 

ct^ °T fr ° m ** ClaSSifiCr "» ** S6C0nd «** the secon 
classifier and determinmg whether the first and second outputs are indicative of an 
acute cardiac ischemic condition. 

15. The device of Claim 13, wherein the cbsaiiier executed by the 
proceasmg unit «*— a « of .oca. features derived from me patient in 
combmauon with the set of global features to determine whether an acute «rdiac 
ischemic condition is detected. 

16. A method of detecting and reporting a condition associated with acute 
cardiac ischemia in a patient, comprising: 

(a) obtaining one or more leads of ECG data from a patient 
ECG data- ^ * ° f ***** fr ° m P atient ' s 

(c) generating a local classification statistic by: 
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(i) deriving a set of local features from the patient; and 

(ii) classifying the set of local features to produce the local 

classification statistic; 

(d) generating a global classification statistic by: 

(i) producing a set of global features by projecting the 
vector of heartbeat data onto one or more basis vectors that define an acute cardiac 
ischemic ECG subspace or a non-ischemic ECG subspace; and 

(ii) .._ classifying the set of global features to produce the 
global classification statistic; 

(e) classifying the local classification statistic and the global 
classification statistic to determine whether the local and global classification statistics 
are indicative of an acute cardiac ischemic condition; and 

(f) reporting whether an acute cardiac ischemic condition is 

detected. 

1 7. The method of Claim 1 6, wherein forming the vector of heartbeat data 
includes: 

(a) analyzing the one or more leads of ECG data to identify one or . 
more heartbeats; 

(b) generating representative heartbeat data for each lead; and 

(c) concatenating the representative heartbeat data for each lead to 
form the vector of heartbeat data. 

18. The method of Claim 16, wherein classifying the set of local features 
includes: 

(a) concatenating the set of local features to form a local feature 

vector, and 

(b) evaluating the local feature vector relative to a representative 
local feature vector derived from a training population to produce the local 
classification statistic. 

19. The method of Claim 16, wherein classifying the set of local features 
includes: 

(a) applying one or more local features in the set of local features 
to a logistic regression equation to produce a probability of detection that is used as a 
composite local feature; 
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(b) dichotomizing the composite local feature; and 

(c) classifying the dichotomized comoosite In^i • 
producing the local classification statistic. ? ^ m 

includef' ^ ° f CIaim ^ ** **■ * of local features 

vector; " *° « ° f W »«« * «» a local feature 

vector and a £ * - — 

classificationstlL ^ * ^ *— h « *. loc, 

21. The method of Claim 20 wherein th* w • , . 

div*ed r a pMty of „ mi . re ^~ ». £ 

local feature vector and te representative ^ "»». 

indudef ^ ° f ^ 16 ' Wttreta the sc. of ,„cal features 

(a) defining a plurality of groups, wherein each group has - 

vector; ^ concatenating the set of local features to form a local feature 

l^t 1 » (C> 8r ° UP ' 3 M^aoobia distance between the 

vector and th. r„uve local feantre vector associated wul Ae 

(d) identifying the local feature vector with one of the nrouns 
ased on fte MahaUnobU distance coated for each group, the group d^dC 
being used as a composite local feature; and P KKntmcauon 

classidcadoniL ^ *** ^ " »"*«*« * 
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23. The method of Claim 16, wherein classifying the set of local features 
includes: 

(a) defining a plurality of groups wherein each group has a logistic 
regression equation associated therewith, the logistic regression equation for each 
group having regression coefficients derived from a logistic regression model using 
selected local features of a training population; 

(b) applying corresponding selected local features derived from the 
patient to the logistic regression equation associated with each group, thereby 
producing for each group a probability of detection that is used as a composite local 
feature; and 

(c) classifying the composite local features in producing the local 
classification statistic. 

. 24. The method of Claim 16, wherein producing the set of global features 
includes: 

(a) calculating an inner product of the vector of heartbeat data and 
one or more basis vectors that define an acute cardiac ischemic ECG subspace to 
produce a corresponding number of ischemic condition projection coefficients; 

(b) . calculating an inner product of the vector of heartbeat data and 
one or more basis vectors that define a non-ischemic ECG subspace to produce a 
corresponding number of non-ischemic condition projection coefficients; and 

(c) using the ischemic condition projection coefficients and the 
non-ischemic condition projection coefficients as the set of global features. 

25. The method of Claim 16, further comprising: 

(a) defining a plurality of groups wherein each group has basis 
vectors associated therewith that define an acute cardiac ischemic ECG subspace and 
a non-ischemic ECG subspace; 

(b) categorizing the patient's ECG data into a group in the plurality 

of groups; and 

(c) using the basis vectors of the group into which the patient's 
ECG data is categorized as the basis vectors onto which the vector of heartbeat data 
is projected. 

26. The method of Claim 25, wherein categorizing the patient's ECG data 
into a group includes: 
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(a) selecting a local feature derived from the patient; and 
selected local feaL. ***** * M ECG d * a *> * « based on the 

27. ThemethodofClaim25 ) furtherco m pri S ing- 
location of JL l^tZ? ^ iS « * 

whet™, category the pata ,'s ECG data too a group include 

* e ST ° f « - - of ft. ,ead s 

calculated; ^ f0min8 " ** f °' Whkh ST *»*- was 

(c) calculating a composite ST elevation for each subgrouo^bv 
^aT" ^ * ST *— - - ^ eacl 

29. The method of Claim 16, further comprising: 

(a) using a Karhunen-Loeve transformation to calculate a first set 
as, vectors *at define an acute cardiac ischemic ECG subspace and a s co d Z 

of basis vectors that define a non-ischemic ECG subspace; 

(b) selecting one or more basis vectors from the first set of basis 
vectors and one or more basis vectors from the second -t «f h ■ 

a L . , . second set of basis vectors as the basis 

vectors onto which the vector of heartbeat data is projected. 

includes 30 ' meth ° d ° f Qaim Whereh daSSif ^ 8 Set of iW-i Matures 

vectored '"""^ 
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(b) evaluating the global feature vector relative to a representative 
global feature vector derived from a training population to produce the global 
classification statistic. 

31. The method of Claim 30, wherein a Gaussian classifier is used to 
evaluate the global feature vector relative to the representative global feature vector 
and produce the global classification statistic. 

32. The method of Claim 16, wherein classifying the local and global 
classification statistics includes: 

(a) producing a combined classification statistic by evaluating the 
local and global classification statistics relative to corresponding representative local 
and global classification statistics derived from a training population; and 

(b) comparing the combined classification statistic with a threshold 
to determine whether an acute cardiac ischemic condition is detected. 

33. The method of Claim 32, wherein the threshold is selected in 
accordance with a desired sensitivity/specificity tradeoff. 

34. A method of detecting and reporting a condition associated with acute 
cardiac ischemia in a patient, comprising: 

(a) obtaining one or more leads of ECG data from a patient; 

(b) forming a vector of heartbeat data derived from the patient ! s 

ECG data; 

(c) deriving a set of local features from the patient; 

(d) producing a set of global features by projecting the vector of 
heartbeat data onto one or more basis vectors that define an acute cardiac ischemic . 
ECG subspace or a non-ischemic ECG subspace; 

(e) jointly classifying the set of global features and the set of local 
features to determine whether the global features and local features are indicative of 
an acute cardiac ischemic condition; and 

(f) reporting whether an acute cardiac ischemic condition is 

detected. 



35. 
includes: 



The method of Claim 34, wherein forming the vector of heartbeat data 
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(a) analyzing the one or more leads of ECG data to 
more heartbeats; a t0 ldentlf y on e or 

00 generatingrepresentativeheartbeatdatafor eachlead- and 
formmeve^ 

(b) dichotomizing the composite local feature. 

to vector denved ton, a traiBng pop^ ^ M . 
used as the composite local feature. 

31 The method of Claim37, wherein the training population has been 

denved for each group, and wherein a Mahalanobis diaaoce is calculated between Z 
tad ta. vector and the represent local feature vector for eTToup 
resuhmg » a plurahtv of composite ,ocal features that ace induded * 
features jointly classified with the set of global features. 

39. The method of Claim 34, wherein the set of local features iointlv 

from, definmg a ptarabty of groups wherein each group has a representative local 

— anng «he set of local features , 0 form a iooaf feature vector, c^g a 
Mafias distance for each group measured between me focal feamre vector 1 
he epresenunv, local feahtre vector associated with ench group, and identifying 
focal feahtre vector with one of the groups hasefi on the lalanobis t£ 
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calculated for each group, the group identification being used as the composite local 
feature. 

40. The method of Claim 34, wherein the set of local features jointly 
classified with the set of global features includes composite local features calculated 
by: 

(a) defining a plurality of groups wherein each group has a logistic 
regression equation associated therewith, the logistic regression equation for each 
group having regression coefficients derived from a logistic regression model using 
selected local features of a training population; and 

(b) applying corresponding selected local features derived from the 
patient to the logistic regression equation associated with each group, thereby 
producing for each group a probability of detection that is used as a composite local 
feature. 

41. The method of Claim 34, wherein producing the set of global features 
includes: 

(a) calculating an inner product of the vector of heartbeat data and 
one or more basis 'vectors that define , an acute cardiac ischemic ECG subspace to 
produce a corresponding number of ischemic conditionprojection coefficients; 

(b) calculating an inner product of the vector of heartbeat data and 
one or more basis vectors that define a non-ischemic ECG subspace to produce a 
corresponding number of non-ischemic condition projection coefficients; and 

(c) using the ischemic condition projection coefficients and the 
non-ischemic condition projection coefficients as the set of global features. 

42. The method of Claim 34, further comprising: 

(a) defining a plurality of groups wherein each group has basis 
vectors associated therewith that define an acute cardiac ischemic ECG subspace and 
a non-ischemic ECG subspace; 

(b) categorizing the patient ! s ECG data into a group in the plurality 

of groups; and 

(c) using the basis vectors of the group into which the patient's 
ECG data is categorized as the basis vectors onto which the vector of heartbeat data 
is projected. 
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into a grip in J" ^ Cldm * Wh6rein ~ S - «» data 

(a) selecting a local feature derived from the patient and 
selected local fl, ^ ^ *» * ^ * ^ b - d « the 

44. The method of Claim 42, former comprising: 

a location of an acute cardiac ischemic condition; and 

W an acute cardiac ischemic condition is detected, reporting the 
locauon of the acute cardiac ischemic condition corresponding to me ™ into 
which the patient's ECG data is categorized. * group into 

45. The method of Claim 44, wherein the location of an acute cardiac 
tschemtc condition is determined based on an ST elevation measured on a lead Id 
wherein categorizing the patient's ECG data into a group includes- 

oh, • At ^ • ST eleVati ° n ° f ° De 0r more of the leads 

obtained from the patient; 

calculated; ^ ^ ** ** ^ ST e,ev ^ WaS 

«1r . (C) u " C ° mP0Site ST devation for each subgroup by 

calculating a mathematical combination of the ST elevation of the leads in each 
subgroup; and 

sub™ u Categ ° rizing *• P atient ' s EC G into a group according to the 

subgroup whose composite ST elevation is greatest. 

46. The method of Claim 34, further comprising: 
of h, • ♦ ^ „ " Karhunen " Loeve transformation to calculate a first set 

a" V TTT " ^ iSCh6miC ECG -i- ** d * -ond set 
of basis vectors that define a non-ischemic ECG subspace; 

(b) siting one or more basis vectors from the first set of basis 
LIT ™^ 

vectors onto which the vector of heartbeat data is projected. 

ft* 4? ', * ^ T ° d ° f Claim 34 ' WhSrein j0iDtly daSS ^ the set of global 
features and the set of local features includes: 
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(a) concatenating the set of local features and the set of global 
features to form a global/local feature vector, and 

(b) producing a global/local classification statistic by evaluating the 
global/local feature vector relative to a representative global/local feature vector 
derived from a training population; and 

(c) comparing the global/local classification statistic with a 
threshold to determine whether an acute cardiac ischemic event is detected. 

48. The method of Claim 47, wherein a Gaussian classifier is used to 
evaluate the global/local feature vector relative to the representative global/local 
feature vector and produce the global/local classification statistic. 

49. The method of Claim 47, wherein the threshold is selected in 
accordance with a desired sensitivity/specificity tradeoff. 

50. A device for detecting and reporting a condition associated with acute 
cardiac ischemia in a patient, the device receiving ECG signals sensed by electrodes 
placed on the patient, the device comprising a processing unit that evaluates ECG 
data acquired from the ECG signals and determines the presence of an acute cardiac 
ischemic condition by analyzing the ECG data, calculating a classification statistic 
reflective of cardiac condition using at least one characteristic obtained from the ECG 
data, and comparing the classification statistic with a threshold, wherein the 
processing unit uses a threshold selected in accordance with a desired 
sensitivity/specificity operating point for the device. 

51. The device of Claim 50, wherein the device further includes a user 
input in communication with the processing unit and wherein the sensitivity/specificity 
operating point of the device is selected by a user of the device via the user input. 

52. The device of Claim 50, wherein the processing unit uses a threshold 
selected for higher specificity for a second evaluation of the patient after the 
processing unit detects an acute cardiac ischemic condition in a first evaluation. 
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